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Abstract. - We study quantum transport of matter waves in anisotropic three-dimensional dis- 
order. First, we show that structured correlations can induce rich effects, such as anisotropic 
suppression of the white-noise limit and inversion of the transport anisotropy. Second, we show 
that the localization threshold (mobility edge) is strongly affected by a disorder-induced shift of 
the energy states, which we calculate. Our work is directly relevant to ultracold-matter waves 
in optical disorder, and implications on recent experiments are discussed. It also offers scope for 
further studies of anisotropy effects in other systems with controlled disorder, where counterparts 
of the discussed effects can be expected. 
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Introduction. — Coherent transport in disordered 
edia is strongly affected by anisotropy effects. It occurs 
in many systems, e.g. electrons in MOSFETs [1], diffusing- 
wave spectroscopy [2], biomedical imaging [3], as well as 
Hght in liquid crystals [4,5], in phosphides [6], or in mi- 
Crocavities [7]. So far, theoretical analyses mainly focused 
on models of disorder made of isotropic impurities imbed- 
ded in anisotropic media [8, 9] or stretched scatterers in 
isotropic media [10, 11], which fairly describe the above 
systems. The recent advent of systems where the disorder 
correlations can be controlled, e.g. tunable arrangements 
of scatterers for microwaves [12,13], engineered optical ma- 
terials [14] or ultracold atoms in optical disorder [15,16], 
opens new perspectives, and it becomes increasingly im- 
portant to better understand transport in disordered me- 
gia with more complex correlations. 

In this letter, we study quantum transport and Ander- 
son localization ( AL) of matter waves in three-dimensional 
(3D) disorder with structured, anisotropic correlations. 
We first determine the diffusion and localization tensors 
using the self-consistent approach of ref . [8] . We show that 
appropriate correlations can lead to rich transport prop- 
erties, such as anisotropic suppression of the white-noise 
limit and inversion of the transport anisotropy. Then, go- 



ing beyond the standard on-shell approximation of ref. [8] , 
wc include the real part of the particle self energy. While 
the latter affects the above results only quantitatively, we 
show that it strongly modifies the behavior of the mobil- 
ity edge. Our results have direct implications to recent ex- 
periments with ultracold atoms in optical disorder [17,18], 
which we discuss. They can also be extended to waves with 
different dispersion relations and other models of disorder, 
where counterpart effects can be expected. 

Quantum transport theory. Let us start with a 
brief theoretical reminder. The building block to describe 
the density (intensity) propagation of a wave in a disor- 
dered medium is the four-point vertex $k k' (q, cj, E') = 

(k+|G(i;+)|k+)(k'_|Gt(£;_)|k_) with G the retarded 
Green operator, k± = k ± q/2 and k'j_ = k' ± q/2 the left 
and right entries, E± = E ± fiaj/2, and (q, to) the Fourier 
conjugates of the space and time variables ^ . The vertex 
$ is governed by the Bethe-Salpeter equation, which can 
be formally written [19] 
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The first term in eq. (1) describes uncorrelated propa- 
gation of the field and its conjugate in the disordered 
medium. The second term involves the vertex function 



-"-Here, we use /(q, uj) = J drdt /(r, t) cxp[— j(q ■ r — ujt)] . 
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U, which includes all irreducible scattering diagrams, and 
accounts for all correlations in the density propagation. In 
the independent scattering (Boltzmann) and weak disor- 
der (Born) approximations [8], \J\^,\^'{q,uj, E) ~ C(k — k'), 
where C (k) is the disorder power spectrum (Fourier trans- 
form of the correlation function (see footnote^)). Here we 
choose the zero of energies such that the disorder is of zero 
average, i.e. V = 0. For any other choice of the energy ref- 
erence all energies appearing below should be shifted by 
V, i.e. replace E hy E —V. At this stage, only the ladder 
diagrams in eq. (1) are retained. It represents an infi- 
nite series of independent scattering events, which leads 
to Drudc-like diffusion. The solution of cq. (1) is then 
dominated in the long time (w — > 0) and large distance 
(q 0) limit by a diffusion pole [19], which reads 



$k,k'(q,w,-B) 



27r S{E-e(k))S{E~e{-k')) 
hNo{E) -iw + q-DB(-E)-q 



(2) 



in the on-shell approximation [such that e(k) = e(k') = E, 
where e(k) is the disorder-free dispersion relation] and 
with No{E) the disorder- free density of states. For mat- 
ter waves in free space (that we consider here), e(k) = 
h?]i^ /2m with m the atom mass. The components of the 
Boltzmann diffusion tensor Db(£') are given by [8] 



1 



2n A 



(3) 



where vi = hki/m is the velocity along axis i, ^ ~ 

;i/27r(C'(fcBk-k'))k'|i5 (withk = k/|k| and = \/2mE/h) 
is the on-shell scattering mean free time, and (...)k|£; = 
/ (2^)'' — e(k)] represents integration over the k- 

space shell defined by e(k) = E. The functions 0^ ^ and 
the real- valued positive numbers A^, indexed by n, are the 
solutions of the integral eigenproblem 



E,k 



= A 
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[8]. In the 



normalized by x(^i5.k<^Ik<^£,k/k|£ 
following, the above equations arc solved numerically to 
determine Db(£') in anisotropic 3D models of disorder. 

Boltzmann diffusion in 3D speckles. Here we 
are particularly interested in the effects of structured cor- 
relations in 3D disorder. In this respect, two configura- 
tions used recently for localizing ultracold atoms are well 
suited. In ref. [17] {single- speckle case), an optical speckle 
disorder is obtained using a single Gaussian laser beam 
of waist w and wavelength Al, propagating along the z 
axis, passed through a ground- glass plate and focused by 
a lens of focal distance / [20]. The disorder correlation 
function C(r) has correlation lengths fT|| = AXi^f'^ /ttw'^ 
along z and a± = Xi^f/mv in the orthogonal plane {x,y) 
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Figure 1: (color online) Disorder power spectrum C(k) for 
the (a) single-speckle and (b) coherent-speckles cases (Fourier 
transforms of the formulas given in the appendix) with the pa- 
rameters of refs. [17,18] {a^^/a± ~ 5.8, and for (b) Al/cx — 
2.16). The functions C(k) are represented as iso-value sur- 
faces (at 2V^a\^) and cuts in the planes defined by the trans- 
port eigenaxes (see text); {ua:,Uy,U2} for (a) and {ux = 
(ua;-Uz)/\/2, uy = Uy,uz = {ux + Uz)/V2} for (b). 



(sec appendix). For 4/ > w, C(r) is elongated along 
z. Then, the disorder power spectrum C(k) is isotropic 
in the {k^jky) plane but significantly shorter along the 
kz axis and forms a plain-wheel shaped structure [see 
fig. 1(a)]. In ref. [18] {coherent-speckles case), the dis- 
order results from the interference of two mutually coher- 
ent and orthogonal speckle fields, propagating along the 
z and X axes, respectively. The power spectrum C(k) 
then shows a complex structure, made of the sum of two 
orthogonally-oriented spectra similar to that of the single- 
speckle case, plus a coherence term (see appendix and 
fig. 1(b)). Most importantly, the latter creates two broad 
structures {bumps), centered on the kx = (k^ — k^)/-\/2 
axis. In the following we call Va = yX7(r = 0) the am- 
plitude and Ea-^ = /m(j\ the correlation energy of the 
disorder. 

Figure 2 shows the components of the Boltzmann dif- 
fusion tensor Db(£') versus the particle energy E, for the 
two above models of disorder. Let us discuss first the 
single-speckle case, which not only forms a useful basis for 
the more complex coherent-speckles case but also shows 
several interesting properties itself. 

Firstly, we find that the scattering time is shorter along 
the symmetry axis z (rj^ ^ < ^ ^) for all values of E. 

This property is simply explained by the fact that C(k) 
has a wider extension in the plane {kx,ky) orthogonal to 
kz than along kz , which offers more scattering channels to 
particles travelling along z. However, as it is well-known, 
the scattering and transport mean free times arc different 
quantities in correlated disorder, due to angle-dependent 
scattering [19,21,22]. In particular, the anisotropy of Db 
can significantly differ from that of Tj^^- Here, we find 
that the anisotropy of Db is reversed with respect to that 
of Tj- p, so that Z?3 > I?b'^, for all values of E'.It is due 



'k,_E' 

to the fact that the orbitals 0" - contributing to in 

eq. (3) are associated to values of X% significantly larger 
than those contributing to D^'^ . 
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Figure 2: (color online) Boltzmann diffusion coefficients along 
the transport eigenaxes (eigencomponents of Db) for the (a) 
single-speckle and (b) coherent-speckles configurations and for 
the parameters of fig. 1. The dotted lines are power-law fits 
(-Dg oc E''" ) to the data in the low and high energy limits. The 
insets show the transport anisotropy factors. 



Secondly, the low-energy behaviour of Db {E) does not 
reduce to a white-noise limit. It is due to the strong 
anisotropic divergence of C(k) in the |k| — !• limit: 
for |k| < cr7\ we have C{k) - £(k)/|k| with £(k) = 



'fc,2)i/2. Using the fact 



we find that ^ and 
is of the form 



exph(a||/2a^)2fc2/(fc2_^fc2)]/(fc2 

that c(k) docs not depend on |k| 

do not depend on E, and that 0^ ^ 

(p"(k)/v^. Then, all terms in eq. (3) are topologically 
unchanged and scale as E. The anisotropy of Db thus 
persists and is constant down to arbitrary low values of 
E, and oc ii^, as observed in fig. 2 for fc^ ^ crj^ 
(i.e. E <C Ecr^)- Conversely, a white- noise limit would 
lead to an isotropic diffusion tensor and Db c>c \/E. In 
other words, in our case, the long-range disorder correla- 
tions lead to the anisotropic suppression of the white-noise 
limit. 

Thirdly, in order to get further insight, it is useful to 
note that, for isotropic systems (see also refs. [23-25]), 
eq. (4) is solved by the spherical harmonics Yj™, and that 
only the first term plus the (p-level) harmonics con- 
tribute to Db in eq. (3). For the anisotropic disorder we 



are considering, we find that the orbitals 



EM 



are topo- 



similar nodal surfaces ^. For <^ ^J^, we find that D^'^ 
is dominated by the first term in eq. (3) and by the Y^- 
likc orbital. For 3> o-J^, the situation changes: while 
D3 is still dominated by the Y]°-like orbital, D^'^ is now 
dominated by the Y^^^-like orbitals with a contribution of 
the Yg^^-like orbitals increasing with E. Then, we find 

_E ^ ^^'^ assuming weak topological change of the 
orbitals and the scaling 1 — cx 1/E (both confirmed 
numerically), we get 0^ oc l/fc^ and D^{E) cx E^^^, as 
observed in fig. 2 for iJ ^ E^^ . This scaling was also 
found for isotropic speckle disorder [23]. Remarkably, in 
spite of the different contributing terms in eq. (3) at low 
and high values of E, the transport anisotropy is nearly 
independent of E [see inset of fig. 2(a)]. 

Let us turn to the coherent-speckles case [fig. 2(b)], 
which is even richer owing to the more complex structure 
of the correlation function. Note first that we recover the 
same general properties as for the single-speckle case, in 
particular the reversed anisotropics of scattering and dif- 
fusion, and the anisotropic suppression of the white-noise 
limit. Here however, the transport eigenaxes are the bi- 
sectors {X, Z} = (x =F z)/-\/2 and the axis Y = y, and 
the anisotropy is significantly smaller. For 2^^ < 3.8crJ^, 
the behavior of Db(£') is governed by the central struc- 
ture of C(k) since, in the on-shell Born approximation, 
a particle of energy E probes C'(k) inside the k-space 
sphere of radius 2kE- Then, the directions X and Z are 
nearly identical but differ from the direction Y and we 
find < ~ Z?f . The most interesting effect appears 
for 2kE > S.SctJ^, i.e. when the bumps of C'(k), located at 
k ~ zbS.ScrJ^kx for the parameters of fig. 1(b), contribute 
to the transport. The scattering time ^ becomes highly 
anisotropic and the orbital dominating is strongly dis- 
torted. As a result, is reduced and the corresponding 
anisotropy factor drops by a factor of ~ 4. This effect 
happens to be strong enough to lead to the inversion of 
the transport anisotropy and we find < < for 
2kE > S.ScrJ^ [sec inset of fig. 2(b)]. 

Localization. — So far, we have discussed incoherent 
Boltzmann diffusion. We now include quantum correc- 
tions due to interference of the diffusing paths by incor- 
porating the maximally-crossed diagrams (Cooperon and 
Hikami terms). It yields the dynamic diffusion tensor 
T>^{uj,E) = Db(£') + AT>{lu,E) with [8] 



AT>{uj,E) 



-Db(S) f dq 



1 



TrhNoiE) J (27r)3 -iw + q.DB(S)-q' 



(5) 



As usual in such approach, the integral in eq. (5) shows an 
ultraviolet divergence, which should be regularized. To do 
so it is usefuU to notice that eq. (5) only depends on Db(£') 
[i.e. in particular, not on C(k)], so that the anisotropy of 
the quantum corrections AD is the same as that of Db [8] , 
and the only characteristic length in the direction i is the 



logically similar to the spherical harmonics, i.e. they show 



■^Note that, in contrast to the isotropic case, here the are not 
dogonoratcd in a given i-like level. 
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transport mean free path 1^{E) = 3y/m/2ED^{E). This 
suggests an eUipsoidal cut-off of radius 1/1^{E) along each 
transport eigenaxis, which we do here It corresponds to 
an isotropic cut-off in the space rescaled according to the 
anisotropy factors of DB(i?). Then, following the stan- 
dard self-consistent theory [26], the above equations are 
solved for D*(aj, E) after replacing DB(i?) by D*(aj, E) in 
the integrand of eq. (5). Hence, as is natural in scaling 
theory and also found with the sigma model for localiza- 
tion [8,27], the anisotropy of D*(a;, E) is the same as that 
of Db(-E) by construction. Proceeding (in 3D) in the long 
time limit (w — >■ 0), a threshold energy Ec appears, so- 
lution of DI^{Ec) = det{DB(i^c)}^/^ = h/Vt^m. For 
E > Ec, D*(ti;,iJ) converges to a real definite positive 
tensor when w — > 0. It describes anisotropic normal dif- 
fusion, characterized by the propagation kernel (2) where 
Db (E) is replaced by the quantum-corrected diffusion ten- 
sor D*(£') = lim^^^o D*('^, -E). For E < Ec, one finds 
D„(uj,E) ^ —iujA{E) for uj ^ 0, where A(£') is a real 
positive definite tensor. It characterizes exponential local- 
ization in the propagation kernel (2), with the anisotropic 
localization tensor hiodE) = ^/A{E). 

Figure 3 shows the components of Lioc (for E < Ec) 
and D* (for E > Ec) for the single-speckle and coherent- 
speckles cases. Due to the preservation of the trans- 
port anisotropy by the quantum corrections and the self- 
consistent procedure, the behavior of Lioc and D* is com- 
pletely determined by that of Db . The anisotropy factors 
of Lioc and D* are nearly independent of E, except for 



Although somewhat arbitrary the factor unity between the cut- 
off radius and 1/1^{E) is justified by the agreement we find with 
another approach in the isotropic case, provided that the real part 
of the self energy is included, as discussed below. 



the inversion of anisotropy of the coherent-speckles case. 
For E < Ec, the anisotropy factors of Lioc(i?) are the 
square roots of those of Db(-E). In the low E limit, we 

find iloc(^) « {D'^/D^)^^^E^^^. When E increases, 
L[^^{E) grows and finally diverges at Ec- For E > Ec, 
the anisotropy factors of Dh<(£') are the same as those 
of DB(i?). The quantum corrections are significant only 
close to Ec. For higher values of E, D*(£') ~ DB(i?), and 
in the high E limit, we have Dl{E) cx (Dl^/ DI^)E^^^ . 

Mobility edge. — The on-shell approach used so far 
is expected to fairly describe the quantum transport prop- 
erties [8,23,26]. It is however limited for quantitative es- 
timates of the mobility edge, which require precise knowl- 
edge of the particle self energy in the disordered medium. 
To estimate the mobility edge, one could in principle 
use the more sophisticated approaches of refs. [24,25,28], 
but the fact that the kind of disorder we are interested 
in is continuous with fine anisotropic structures makes 
them hardly practicable numerically. In order to over- 
come this issue, we propose an alternative method based 
on the assumption that the leading term missing in the 
on-shell approximation is the real part of the self energy, 
T,'{k,E) = Pjj0y, ^E^-Z^? (with P the Cauchy princi- 
pal value), which produces a shift of the energy states. We 
incorporate S' (k, E) into the theory by averaging, in first 
approximation, its k-angle dependence. It amounts to re- 
place the on-shell prescription by e(k) = E' = E — A(i?) 
with 

A{E)^^f dn^^'{k,E) (6) 

where is the k-space solid angle. Within this approach, 
all previous quantities are now regarded as functions of 
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Figure 4: (color online) On-shell {E'^) and corrected (E^) mobil- 
ity edges versus the disorder amplitude Vn for the single-speckle 
and coherent-speckles cases and the parameters of fig. 1. 

E' instead of E. In particular, the mobility edge Ec is 
the solution of Ec — A{Ec) = E'^, where E'^ is determined 
using the on-shell approach. Here two remarks are in or- 
der. First, when applied to isotropic disorder, our method 
gives good agreement (within 5 — 7%) with the values of 
Ec found in ref. [25] by a different approach, which in par- 
ticular does not require the introduction of an ultraviolet 
cut-off. Second, the k-angle averaging used in eq. (6) is 
justified a posteriori by the weak k-angle variations of E' 
around its average value at the mobility edge Ec (with 
standard deviations less than 10 — 15%). 

For the single-speckle and coherent-speckles cases, the 
on-shell (E'c) and corrected {Ec) mobility edges are shown 
in fig. 4. It is eye-catching that the shift of the energy 
states completely changes the behavior of the mobility 
edge. While E'c is positive (i.e. above V) and increases 
with Vr, we find that Ec is negative (i.e. below V) and 
decreases with V^. For < Ea^, this behavior is qual- 
itatively similar to that found for isotropic disorder [25]. 
For larger values of V^, Ec further decreases, consistently 
with the expected behavior that Ec should approach the 
percolation threshold deep in the classical disorder regime 
(Fr»S,^) [29]. 

Discussion. — By considering matter waves in two 
particular models of optical disorder, we have shown that 
structured disorder correlations in 3D can induce rich 
quantum transport and AL properties, e.g. anisotropic 
suppression of the white- noise limit, and inversion of the 
transport anisotropy with energy. Such effects can be 
expected for other kinds of waves and/or other models 
of disorder, and are particularly relevant to new systems 
where the disorder correlations can be controlled [12-18]. 
In the case of ultracold atoms, to which our study di- 
rectly applies, the transport properties can be probed 
by direct imaging of the atoms and control of the en- 
ergy distribution. This control does not need to be very 
fine. For instance, in the single-speckle case, we found 



almost constant anisotropy factors {D^/D^'^ ~ 10 and 
^foc/^k)c — •^•2) I Slid experimental data can be com- 
pared to these predictions almost independently of the 
energy distribution. Although no precise value has been 
extracted from the experiment of ref. [17], the experimen- 
tal data indicate significantly larger anisotropy. Further 
analysis would be required to clarify the origin of such 
a discrepancy. Conversely, in the coherent-speckles case, 
the anisotropy factor was shown to be in fair quantitative 
agreement with theory [18]. The inversion of the trans- 
port anisotropy was however not observed because the im- 
ages were taken in the (y, z) plane. It only gave access to 
Dy = and D'- = {D^ + D^)/2, which do not show the 
inversion. In order to observe it, it is required to image 
the atoms along the transport eigenaxes and to tune the 
balance between the populations of low- and high-energy 
states. 

Another major challenge is the evaluation of the 
mobility edge in 3D disorder with structured correla- 
tions. In this case numerical evaluations in advanced ap- 
proaches [25, 28] are hardly practicable. Here we have 
proposed and used an applicable approach. In particular, 
it yields a mobility edge that is negative (i.e. below V) , as 
also predicted for isotropic speckle potentials [25]. Com- 
paring to ref. [17], our calculations hence significantly dif- 
fer from experimental values (e.g. for Vr = 600 nK x feg — 
7.1^;^^, we find Ec ~ -300 nK x while -1-900 nK x 
is measured). However, the method used in ref. [17] to 
infer Ec from the localized fraction uses the free-space 
kinetic energy distribution, which neglects the disorder- 
induced distortion of the energy distribution. It is ques- 
tionable because the latter is, in particular, necessary to 
account for negative energy states (i.e. below the disor- 
der average value). Comparing to ref. [18], we find that 
A (E'c) as calculated here is of the same order of magni- 
tude as the heuristic shift introduced in ref. [18] (e.g. for 
Vn = hx 680 Hz ~ 0.35£:^^, we find A{Ec)/h = -390 Hz 
and the heuristic shift is —225 Hz). Given systematic un- 
certainties in the analysis of the experimental data, it is a 
fair agreement. A precise test of the present theory would 
require a reliable determination of the energy distribution 
in ultracold-atom experiments, which is not available so 
far. 

Beyond these direct applications, our results pave the 
way to further studies of anisotropy effects in coherent 
transport and AL. On the one hand, even more com- 
plex correlations can be designed in ultracold-atom ex- 
periments in 3D [30] as well as in 2D [31,32]. On the 
other hand, it would be interesting to explore counter- 
parts of the discussed effects for waves with different dis- 
persion relations and/or in other kinds of controlled dis- 
order [12-14]. /,From a theoretical point of view, several 
extensions of our approach may be envisioned. On the 
one hand, it would be worth going beyond the shifted on- 
shell approach used here, and incorporate the full struc- 
ture of the particle spectral function. On the other hand, 
it would be interesting to estimate possible corrections to 
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the anisotropy factors in the locahzed regime, which may 
be expected to be significant in strongly anisotropic dis- 
order. This would require developing an approach where 
the preservation of the Boltzmann transport anisotropy in 
the quantum corrected tensor is not built in the theory, as 
in the self-consistent approach used here. 
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APPENDIX 

Correlation functions. — In order to calculate 
the autocorrelation function of the disorder, C(r) = 
V{r + r')V{r'), for the single-speckle and coherent- 
speckles cases, we use the paraxial approximation, which 
yields explicit formulas. For the single-speckle case, we 
find C(r) = V^cisp{x,y,z) with 



cisp{x,y,z) 



1 



1 + 4z2/crj 



■ cxp 



1 + 4z2/ct2 



For the coherent-speckles case, we find C(r) = (Vr/2)2 x 
{cisp(a^, y, z) + cisp(z, y, x) + 2ccoh(a;, y, z)} with 



Ccoh(r) ^J cisp{x, y, z) x cis-p{z, y, x) 



(l + 4f|)cos[0(r)]-l-2^ 



1 -H4zVcr2 /l-|-4x2/cr 



and 



2lT 

Al 



{x 



X +y 
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